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1 . INTRODUCTION 

Many problems in geophysics involve different areas of investigation 
which are interrelated. One such problem is investigation of phenomena at 
tidal frequency in the solid earth and ocean. On the one hand, ocean 
generated tidal signals are present in many geophysical measurements, such 
as those of satellite elevation, gravity, tilt, and strain. On the 
other hand, the ocean tide is affected by the elastic properties of the 
solid earth, and so accurate modeling of the ocean tide depends on knowledge 
of the solid earth. 

In order to make full use of geophysical data, geophysicists and geod- 
esists need realistic models of the open ocean tide. But knowledge of the 
ocean tide has h'^en difficult to obtain because ocean basins are close to 
resonance at tidal frequencies. This makes accurate modeling of the ocean 
tide difficult. It is the purpose of this paper to present models of the 
M2, S2, and K1 tides that were developed to provide a first order correction 
to the problem of resonance. To facilitate the use of these models by 
geophysicists, maps of the geocentric tide, the induced free space potential, 
the induced vertical component of the solid earth tide, and the induced 
vertical component of the gravitation field are presented. These fields 
were calculated with Green's functions provided by W. E. Farrell (private 
communication) . 

2. THE PROBLEM 

Laplace's Tidal Equations (LTE) may be written in negative mercator 
coordinates as 
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- (2n tanht)v = -g(?-r/g)^/a 

H- (2n tanht)u = -g(?-r/g)^/a 

Cc^-6^) H- ({uD)^ + (vD)^)/a sech^T = 0 

where {-x,i>) are the mercator latitude and longitude, (u,v) are the corres- 
ponding velocities, ^ is the geocentric ocean tide, r is the total tide 
generating potential, 6 is the geocentric solid earth tide, D is the local 
depth of the ocean, a is the radius of the earth, and n is the earth's 
angular rate of rotation, t = -an tan(Tr/4-0/2) where 0 is latitude. 
Introducing the observed tide Cq = C-5 and writing 

CoU)T,t) = Re{CQ(^.T)e“^''^} etc. 

allows LIE to be written as 

-iau - (2n tanh.T}v = -g(tQ+<s-r/g)_j_/a 

-iav + C2n tanht)u = -g(Cp+<s-r/g)ya 

-icr?Q + + (vD)^)/a sech^t = 0. 

Solving for u and v, and placing these results back into the continuity 
equation yields 

jCCCq) + sech^T =;£(r/g-6) 

where 

^ = QHv^ + [(QH)^ - Ci/s)(QH tanhx)^] 3^ 

+ [CQH)„ + (i/s)(QH tanhT}^]3^ 
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Q = l/(s^ “ tanh^t) 


= 4f}2a2/gD^ 

D = D^H(«jt,T) s = o/2n. 

For a rigid earth, r = U (the astronomical potential), and 6=0. 

In the presence of solid earth deformation and self-gravitation by the 
water column, r/g-6 may be expressed after Hendershott (1972) as 

r/g-6 (l+kp-hj^)U^/g ,t')G(<}»,t/(|i' ,T') sech^T'dT'd<|t' , 

where 6(4>,T/if>' , t‘ ) is a Green's function constructed by Farrell (1972), 
(kn»hj^) are appropriate Love numbers (Munk and McDonald, 1960), and is 
the n-th spherical harmonic of U. 

Thus the final equation becomes 

oce*ans 

sech^T =JC[XO+k^-h )U^/gl , 
n 

We call this the elevation equation. Regardless of boundary conditions and 

regardless of the (numerical) procedure used to invert the operator.^, 
the presence of the global integral in the above equation makes solution 

difficult. 

All this paper's models specify the observed tide at the computational 
boundaries, as in Hendershott (1972). The models are frictionless but 
energy dissipation is obtained by allowing energy to flow through the 
computational boundaries, thus modeling energy flow into shallow marginal 
seas and shelves. 
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Hendershott (1972) proposed the iterative sequence 


+ G^sech^ =X[S d+k„-h^)U^/g] 

n 

^ ^ ^Ssech^T'dT'd<{i' J 


oceans 


and found it to be divergent for M2 in the absence of interior dissipation. 
Gordeev, ^ aj^. (1977) showed that this procedure will converge in the 
presence of interior dissipation. 

Parke (1978) used iterates generated as above to provide a basis 

set for a least squares solution to the dissipationless equations. The 
iterative procedure was restarted periodically to avoid redundancy in the 

These least square solutions, however, proved unrealistically resonant 
for the semi-diurnal constituents M2 and S2. None of the solutions conserved 
mass, even that for K1 (which compared favorably with island data). This is 
an important liability for geophysical calculations, as was shown by Farrell 
(1972a). 


The difficulty with being near resonance is that small errors in how 
one's model represents the real ocean cause small errors in the frequencies 
at which the model basin resonates. Far from resonance, these errors should 
have a correspondingly small effect on the model tide. Near-resonance, 
however, these small errors cause significant errors in the amplitudes of 
near resonant modes, and consequently cause significant variations in the 
model tide. Even relatively small changes in basin parameters, such as the 
mean depth of the ocean, cause significant changes in the represented tide 
near resonance. 
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In a one mode system, i.e., one where a single near-resonant mode 
dominates the solution, the amplitude of the mode can be adjusted by varying 
a single parameter such as the mean depth. A good example of this is the 
model of the Gulf of California by Stock (1976), where empirically choosing 
this single parameter improves the model dramatically,’ Unfortunately, the 
global tide appears to be a multi-mode system, and so adjustment of the tidal 
elevation field via a parameter search would require at least two and most 
likely more parameters. As the cost of such a search depends quite strongly 
on the number of parameters used, this approach could prove quite expensive. 

An alternate method for adjusting the amplitude of near resonant modes 
is given in the next section. 


3. A DYNAMIC INTERPOLATION OF ISLAND DATA 

The test functions used by Parke (1978) to solve LTE in the 

presence of tidal loading may be regarded as linearly independent combina- 
tions of the eigenmodes of the discretized tidal operators,. We have 
generated a representation of the actual ocean tide by using these test 
functions to interpolate between island data in the least squares sense. 

If we believe that the modifications of our finite difference operator 
which are necessary to make its near tidal resonant frequencies the same as 
those of the real ocean would have but small effect on the spatial shape of 
the corresponding eigenmodes, then this least squares fit closely approximates 
the solution of that modified finite difference model. 


A 

The least squares fit 


n 

1 A. 
i=l ^ 



is 


obtained by solving 
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(9/8 A, ) ( S lie,- X A, ;('>(♦[, T,)]|^) = 0, J=l.n 
^ r=islands ' 1=1 ' “ ‘ ‘ 

for the weights Here Cj is the tide at island I, 

The test functions each separately reproduce the specified 

n 

boundary tides. Requiring S A. =1 forces c to also reproduce them. The 

\j ^ 

test functions do not individually conserve mass, i.e., 

^ 0 ^^^ - ^ Sf sech^Tdtd^l/iarea of oceans}?^ 0. 

oceans 

Requiring XAj = 0 forces r,^ to conserve mass. 

For the M2 fit, n=24 iterates were used, for S2, n=16 and for K1 , n=8. 

n n 

The fits presented below were obtained requiring ^ {although SAj 
came out within a few percent of unity when this constraint was not imposed). 
^“0 was likewise required although it also tended to be small (less 

J C 

than .1 cm for all three constituents when unconstrained). We do not entirely 

A 

understand why thus conserves mass so well while the least squares solutions 

to the finite difference elevation equation, with the same boundary conditions, 
n \ 

do not (2 Aj ^ was several orders of magnitude larger for the least 
square solutions even when the mean depth was intentionally altereo simply 
to avoid resonance). 

Least squares fits using different sets of M2 test functions generated 
over a range of 5% in the mean depth yield essentially the same elevation 
fields. This occurs even though the different sets of test functions differ 
markedly. We regard this latter point as qualitative confirmation of our 
belief, expressed above, that changes in the spatial shapes of the eigenmodes 
of the tidal operator are not a major source of error in our calculations. 
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Island tide gauges tend to be in the worst possible places, i.e., 
inside lagoons. This problem can cause a bias in the generated elevation 
field. To get an idea of this bias, seven islands in the Pacific were 
empirically corrected for near station phrise lag by means of tsunami travel 
time charts, as was suggested by Mr. B, D. Zetler. The corrections for M2 
are as follows; Canton Island 0 , Apra K:,rbor (Guam) 5 , Johnston Island 
7 , Kwajalein 5 , Midway Island 4 , Tahiti 5 , Honolulu 5 . These phases 
are to be subtracted from the published M2 phases for these stations. 
Corrections for S2 and K1 are of similar small magnitude, 

Larsen (1977) has shown how to remove the effect of extended island 
arcs (i.e., the Hawaiian Chain) from ir^land tide observations. His method 
should be applied systematically in future refinements of the present study. 

It should be noted that there is a growing set of coastal and deep sea 
pressure gauge measurements, which would provide an ideal data set for 
these calculations. At the present time, these measurements are being 
collected for publication under the auspices of the International Association 
for the Physical Sciences of the Ocean (lAPSO) by Cartwright and Zetler 
(Zetler, personal communication). 

The M2 fit to observed tides for 3930 m mean depth plus scatter plots 
comparing It vvith island data are given in Figures 1 and 2. The least 
squares M2 elevation field shows structure similar to that of the solution 
of Zahel (1977). The main deviation occurs in the South Atlantic where 
Zahel shows an amphidrome not found in the least squares field. Otherwise, 
Zahel 's amplitudes in the Pacific are significantly larger and his phases 
are entirely different there. The elevation field in the Atlantic and 
Indian Oceans is likewise similar tc the corresponding cotidal charts of 
Dietrich (1944). 
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The S2 observed elevation field and corresponding scatter plots are 
given in Figures 3 and 4. It should be noted that the S2 tide includes 
significant radiational as well as gravitational forcing. Zetler (1971) 
indicates that the amplitude ratio between the radiational and gravitational 
tides off the coasts of the US is 16%. Since both the coastal and island 
elevations reflect this dual forcing, the least squares solution is expected 
to fit island data well even though the test functions were generated with 
only gravitational forcing. 

The least squares K1 observed elevation field is similar to that of 
Dietrich (1944) throughout the world's oceans. The elevation field and 
its associated scatter plot are given in Figures 5 and 6. 

Table 1 summarizes some features of the models, including the rate of 
working by the tide generating forces, the global average potential energy, 
and the global Q, all estimated as in Hendershott (1972); kinetic energy 
calculations involved unresolved numerical problems due to the proximity of 
the inertial latitudes to mesh points and hence equipartition was assumed. 
Platzman's (1975) results suggest that this involves an error of only a 
few percent. 

19 

The calculated rate of lunar working for the M2 model, 2.22 x 10 
1 Tg 

erg s" , is smaller than the usual estimates, which range from 3.04 x 10 
erg s~^ for the Hendershott (1972) model to 3.77 x 10^^ erg s ^ for Zahel 
(1977). This smaller value is in better agreement with the estimate of the 
M2 dissipation in shallow seas by Miller (1966) of 1.7 x 10 erg s . 

Astronomically generated estimates currently seem to be in a state of flux. 

19 -1 

Estimates recently have risen from the long accepted 2.7 x 10 erg s 
(derived from the Spencer~Jones value for tne deceleration of the lunar 
longitude) and have subsequently fallen again. Muller (1976) reviews them 
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19 -1 

and places the present value at 3.3 + .2 x 10 erg s . Goad and Douglas 
(1978) obtain a very similar value fror analysis of satellite orbital param- 
eters. It should be noted that this estimate represents the sum of all lunar 
dissipationso and as such should be considered simply an upper boundary for 
the M2 tide. Estimates of the non-M2 oce^n effect on the lunar dissipation 
vary from 0.3-1 ,0 x 10^^ erg s“\ thus reducing the above figure to 2. 1-3.2 x 
10^ ' erg s"^ . 

The nearness to resonance of the semi-diurnal constituents M2 and 32 
makes it difficult to model these tides. The global Q values of 17.0 for 
M2 and 30.1 for S2 are consistent with the nearness to resonance. The K1 
constituent on the other hand has every appearance of being far from reson- 
ance. The basic iterative procedure of Hendershott (1972) converges, and 
there is little difference between this solution and the least squares fit 
to island data. The only significant advantage of the fit to island data 
is that it better conserves mass. The apparent distance from resonance 
occurs despite the fact that PI atzman (1975) finds several resonances near 
the diurnal frequency. No satisfactory explanation has been found. The 
calculated Q value for the island fit representation of 5.8 is consistent 
with the paradoxical absence of resonance. 

4. MAPS OF THE GEOCENTRIC OCEAN TIDE, THE OCEAN INDUCED POTENTIAL, THE 
OCEAN INDUCED VERTICAL COMPONENT OF THE SOLID EARTH TIDE, AND THE 
OCEAN INDUCED GRAVITY TIDE 

These maps are presented for their geophysical interest. The geocentric 
ocean tide ? = + 6 is directly observable by satellite altimetry. The 

total tide generating potential at the mean sea surface -is 
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r = (HI< 2 )U 2 + S(l+k'^)ga^ Cgp. 


the second half of this expression is the ocean induced potential at the 

mean sea surface. When analytically continued upward, kgUg + 2(l+k'|j)ga^ 

perturbs satellite orbits. It is not plotted since It may be obtained by 

/ 

adding k 2 U 2 to the previous plot. The vertical solid earth tide is 
5 = h 2 U 2 /g + 2 h'^ the second half of this expression is the ocean 

Induced vertical component of the solid earth tide. The ocean induced gravity 
tide as measured by a gravimeter on the solid earth's (moving) surface is 

g/aS (n+2 h'„ -Cn+Dk'^) 

Here h'^^ and k'^ are loading Love numbers (Munk and McDonald, 1960), 

= (3/2n+l ) (water density/earth density), is the n-th spherical 
harmonic of the observed ocean tide and a is the radius of the earth. 

Note that the ocean induced potential at the solid earth's moving 
surface (i.e., not at mean sea level) is 2(l+k'^-h'^)ga^ It is 

readily obtained by differencing the ocean induced potential at mean sea 
level and g times the ocean induced vertical component of the solid earth 
tide. It is therefore not plotted. 

The maps are pr'^sented in Figures 7 tlirough 18 in the order given in 
the title of this section. For each physical quantity the maps are in the 
order M2 then S2 then Kl. Numerical tables of the mapped quantities are 
presented in a microfiche appendix available in limited supply from the 
office of Marine Geodesy, 
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It should be noted that these calculations do not include the ocean 
tides In coastal seas and shelves. Thus if one is interested in gravity or 
solid earth tides in a region near the coast, a local correction must be 
ii'ade using observations or reconstructions of nearby offshore tides. 

All calculations were performed with Green's functioris calculated by 
Farrell (1972). 

Results of the K1 calculation are of particular interest because there 
is a proposed core resonance of near-Kl frequency (Rochester, 1973). Maps 
of the induced component of the gravitational field can therefore be used 
to indicate regions in which there is a likelihood that the Kl ocean 
contribution will be small; these then are optimum locations at which to 
look for the core resonance in gravity measurements. Indications are, from 
the calculations made here (Figure 18), that the best place for such measure- 
ments would be in Northern Russia where there is a large region with very 
small amplitudes evidently due to the great distance from the oceans. 

Because of the great extent of this region, positional errors should have 
little Impact. The next most likely place would be in Australia where there 
appears to be an amphidrome. Unfortunately, the gradient away from the 
amphidrome is large, making cor»'ect placement of the amphidrome critical to 
finding the low amplitude region. In addition there is apparently an 
amphidrome in Antarctica, but as the core resonance signal dies off toward 
the pole, this may not be an optimal region for its measurement. 
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FIGURE CAPTIONS 


Figure 

Figurii 

Figure 
Figure 
Figure 
Figure 
Fi gure 

Figure 
Figure 
Fi gure 

Figure 

Figure 

Figure 


1. Dynamic interpolation of the global M2 tide based on test 
functions generated at 3930 m mean depth (amp. in cm, phases in 
degrees relative to passage of the tide generating body over 
Greenwich). 

2. Scatter diagrams comparing modeled M2 tidal elevations with 
island data. 

3. The same as Figure 1, except for S2. 

4. The same as Figure 2, except for S2. 

5. The same as Figure 1, except for Kl. 

6. The same as Figure 2, except for Kl . 

7. The geocentric M2 tide based on the model of Figure 1 (units 
as in Figure 1). 

8. The same as Figure 7, except based on the S2 model of Figure 3. 

9. The same as Figure 7, except based on the Kl model of Figure 5. 

10. The free space potential (divided by g) induced by the M2 model 
of Figure 1 (amp. in cm, phase in degrees as in Figure 1). 

11. The same as Figure 10, except induced by the S2 model of Figure 3. 

12. The same as Figure 10, except induced by the Kl model of Figure 5. 

13. The vertical component of the solid earth tide induced by the 

M2 model of Figure 1 (amp. in cm, phases in degrees as in Figure 1). 
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Figure 14. 
Figure 15. 
Figure 16. 

Figure 17. 
Figure 18. 


The same as Figure 13, except induced by the S2 model of Figure 3. 

The same as Figure 13, except induced by the K1 model of Figure 5. 

The vertical component of the gravitational field induced by the 
M2 model of Figure 1 (amp. in 10 cm s , phases in degrees 
as in Figure 1 ) . 

The same as Figure 16, except induced by the S2 model of Figure 3. 

The same as Figure 16, except induced by the K1 model of Figure 5. 
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